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ABSTRACT 

We compile a holistic scenario for molecular cloud (MC) evolution and control of the star formation 
efficiency (SEE), and present a first set of numerical tests of it. A lossy compressible cascade can 
generate density fluctuations and further turbulence at small scales from large-scale motions, implying 
that the turbulence in MCs may originate from the compressions that form them. Below a sonic scale As, 
turbulence cannot induce any further subfragmentation, nor be a dominant support agent against gravity. 
Since progressively smaller density peaks contain progressively smaller fractions of the mass, we expect 
the SFE to decrease with decreasing As, at least when the cloud is globally supported by turbulence. 
Our numerical experiments confirm this prediction. We also find that the collapsed mass fraction in 
the simulations always saturates below 100% efficiency. This may be due to the decreased mean density 
of the leftover interclump medium, which in real clouds (not confined to a box) should then be more 
easily dispersed, marking the "death" of the cloud. We identify two different functional dependences 
("modes") of the SFE on As, which roughly correspond to globally supported and unsupported cases. 
Globally supported runs with most of the turbulent energy at the largest scales have similar SFEs 
to those of unsupported runs, providing numerical evidence of the dual role of turbulence, whereby 
turbulence, besides providing support, induces collapse at smaller scales through its large-scale modes. 
We tentatively suggest that these modes may correspond to the clustered and isolated modes of star 
formation, although here they are seen to form part of a continuum rather than being separate modes. 
Finally, we compare with previous proposals that the relevant parameter is the energy injection scale. 



1. INTRODUCTION 

Two of the main questions about molecular cloud (MC) 
structure and star formation concern the origin of molec- 
ular cloud turbulence and of the low efficiency of star for- 
mation (SFE). Indeed, the turbulent properties of molecu- 
lar clouds are remarkably uniform, even in clouds without 
stellar energy sources (see, e.g., McKee 1999 and refer- 
ences therein). And, defined as the fraction of the mass in 
a MC that ends up in stars before the cloud is dispersed, 
the SFE is well known to be as low as a few percent with 
respect to the total cloud mass (see, e.g., Evans 1991 and 
references therein). 

It is known that turbulence can provide additional sup- 
port against global gravitational collapse (Chandrasekhar 
1951; Bonazzola et al. 1987, 1992; Vazquez-Semadeni & 
Gazol 1995). In particular, numerical studies have sug- 
gested that turbulence is able to prevent the collapse 
when the energy injection scale is smaller than the Jeans 
length at typical core densities (Leorat, Passot & Pou- 
quet 1990; Klessen, Heitsch & Mac Low 2000, hereafter 
KHMOO). Simultaneously, supersonic, compressible turbu- 
lence is known to produce "fragmentation" , i.e., the forma- 
tion of smaller-scale density substructures that can them- 
selves possibly undergo local collapse (Sasao 1973; Hunter 
& Fleck 1982; Tohhne, Bodenhcimer & Christodolou 1987; 
Scalo 1987; Elmegreen 1993; Padoan 1995; Vazquez- 
Semadeni, Passot & Pouquet 1996; KHMOO; Padoan et 
al. 2001; Heitsch, Mac Low & Klessen 2001). Padoan 
(1995) and KHMOO have proposed that turbulent fragmen- 
tation has an important role in determining the SFE. The 



latter authors have shown that the efficiency in numeri- 
cal simulations (measured as the fraction of mass in col- 
lapsed objects) increases as the energy injection scale is in- 
creased. However, this result was obtained in simulations 
with an idealized random forcing applied at specific scales, 
and ubiquitously in time and space. Instead, the turbu- 
lence in real MCs is probably injected over a wide range of 
scales, but at localized sites and times (Scalo 1987; Nor- 
man & Ferrara 1996; Avila-Reese & Vazquez-Semadeni 
2001). In this case, the driving scale is not uniquely de- 
fined. Also, if molecular clouds form through the col- 
lisions of larger-scale streams in the more diffuse ISM 
(Elmegreen 1993; Ballesteros-Paredes, Vazquez-Semadeni 
& Scalo 1999; Vazquez-Semadeni, Ballesteros-Paredes & 
Klessen 2002, hereafter Paper I), then they can be part of a 
turbulent cascade originating from much larger scales, and 
thus the driving scale can be much larger than the Jeans 
length in the molecular clouds (Mac Low & Ossenkopf 
(2001). This appears inconsistent with the requirement 
of a driving scale smaller than the cloud's Jeans length in 
order to maintain a low efficiency. 

In Paper I, we have compiled a holistic scenario in which 
both the turbulence in molecular clouds and the SFE 
are part of the life cycle of MCs: MC turbulence origi- 
nates from the very process that forms the clouds, namely 
colliding streams in the diffuse ISM (Elmegreen 1993; 
Ballesteros-Paredes et al. 1999), through bending- mode 
instabilities (Hunter et al. 1986; Vishniac 1994; Walder & 
Folini 1998; Klein & Woods 1998). In this sense, MC tur- 
bulence is part of a compressible and lossy cascade which 
may originate at intermediate scales in the Galactic disk 
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(> 1 kpc) from various instabilities (Parker 1966; Toomre 
1964; Elmcgrcen 1991), supernova energy input (e.g., Mc- 
Cray & Snow 1979; MacLow 2002), and stellar winds (Nor- 
man fc Silk 1980), and continues towards smaller scales, 
as the internal turbulence produces yet smaller-scale tur- 
bulent compressed layers (Kornreich & Scalo 2000). The 
compressible cascade has to be lossy because shocks are 
known to be dissipative by transferring a fraction of the en- 
ergy involved directly to the dissipative scales (Kadomtsev 
& Petviaslivili 1973). This compressible-cascade regime, 
however, must end at a "sonic" scale, denoted Ag, below 
which the turbulence internal to the new structures be- 
comes subsonic and therefore cannot produce any further 
subfragmentation (Padoan 1995). This docs not neces- 
sarily imply that the dissipative (or "inner" ) scale of the 
turbulence has been reached, but only that at this point 
the cascade can turn into an incompressible one. Thus, in 
this scenario, subsonic "cores" arc the natural end of the 
compressible part of the cascade in the ISM. 

When the turbulence becomes subsonic;, it ceases to be 
a stronger source of global support than thermal pres- 
sure, and the Jeans criterion can be applied to deter- 
mine whether a "core" is gravitationally unstable and 
should proceed to collapse (Padoan 1995), or else re- 
bound and merge back into its environment (Elmegreen 
1993; Taylor, Morata & Williams 1996; Ballesteros- 
Paredes et al. 1999; Vazquez-Semadeni, Shadmehri & 
Ballesteros-Parcdes 2002). This effectively sets an upper 
limit to the SFE, because, if globally the cloud is sup- 
ported by its internal turbulence, only massive enough 
clumps among those with sizes / < Ag are susceptible 
to local collapse (Padoan 1995). In this sense, star- 
forming cores constitute the intersection of the set of sub- 
sonic regions with the high-mass, super-Jeans tail of the 
mass distribution (Padoan 1995; Ballesteros-Paredes & 
Vazquez-Semadeni 1997; Shadmehri, Vazqucz-Scmadcni 
& Ballcstcros-Paredes 2001). In particular, Padoan (1995) 
gave a scmi-phenomenological theory for the efScicncy and 
mass distribution of proto-stellar cores in turbulent molec- 
ular clouds (recently extended by Padoan & Nordlund 
2002). In that work, however, he considered only a power- 
law turbulent spectrum, neglecting the necessary existence 
of inner and outer scales, and the role of the energy injec- 
tion scale, which was proposed to be crucial by Leorat et 
al. (1990) and KHMOO. 

In this Letter, besidc^s putting together and advancing 
the above global scenario, we present (§2) a first set of 
numerical tests supporting it, by showing that the scale 
Ag indeed plays a crucial role in limiting the SFE. Clearly, 
this is only a partial ingredient, as the mass distribution 
in the flow is also expected to be crucial (Padoan 1995). 
Our experiments also support the notion of the dual role 
of turbulence in providing support to large scales while 
promoting small-scale collapse, as we discuss in §3, where 
we also give a re-interpretation of the results of KHMOO 
in the framework of our scenario. 

2. NUMERICAL EXPERIMENTS AND RESULTS 

We consider a set of 12 numerical experiments similar 
to those reported by KHMOO. We refer the reader to that 
paper for details (see also Klessen &: Burkert 2000, 2001). 
Here we just mention that they are smoothed particle hy- 



drodynamics (SPH) simulations of non-magnetic turbu- 
lent flows using 205,379 particles, with self-gravity, and 
random forcing applied only at wavenumbers in a narrow 
interval fed — 1 < |k| < fed, where fed = L/ Ad is the driving 
wavcmimbcr, L is the computational box size and Ad is 
the driving scale. The simulations are carried out in non- 
dimensional units, in which there are 64 thermal Jeans 
masses (with respect to the mean density) in the compu- 
tational box, the time unit is to = ts/1.5, where tg is the 
free fall time, the sound speed is c = 0.1, and the mean 
density is po = 1/8. A constant kinetic energy input rate 
is used, with an amplitude chosen to maintain a fixed level 
of turbulent kinetic energy Et. For our analysis, the SPH 
data are interpolated onto a cubic grid of 128^ cells, but 
note that the effective resolution of the SPH method in 
the density peaks (cores) is higher than that. 

The runs are labeled mnemonically as MxKy, where x 
denotes the rms Mach number and y denotes the forcing 
wavenumber fed. The rms Mach number Ms and the energy 
injection scale Ad were taken by KHMOO as the controlling 
parameters in their investigation of the amount of mass 
in collapsed objects M* as a function of time for several 
simulations. The numerical scheme allows this mass to be 
measured through the usage of "sink" particles, which are 
particles that replace a collapsing gas clump once it be- 
comes too dense and small to be properly resolved. A sink 
particle inherits the combined masses of the replaced SPH 
particles, as well as their linear and angular momenta, and 
has the ability to continue accreting further SPH particles. 

In order to investigate the role of the sonic scale Ag, 
we consider the velocity dispersion-vs.-size relation, to- 
gether with the mass-accretion histories of the simula- 
tions. The velocity dispersion-size relation we consider 
is not restricted to the densest, or highest- column den- 
sity projected structures, as is customary in both oser- 
vational (Larson 1981; see Blitz 1993 for a more recent 
review) and numerical (Vazquez-Semadeni, Ballesteros- 
Paredes & Rodriguez 1997; Ostriker, Stone & Gammie 
2001; Ballesteros-Paredes & Mac Low 2002) astrophysical 
studies. Instead, we consider a subdivision of the whole 
computational box in regions of size Z, and measure the 
velocity dispersion in all such regions. The mass accretion 
histories are computed as the fraction of the total mass in 
collapsed objects (sink particles) as a function of time. 

Figure 1 shows the velocity dispersion-size relation for 
all runs. For every run, the solid line shows the average 
velocity dispersion for all regions of size I; the "error" bars 
indicate the range of values of the velocity dispersion ex- 
hibited by all regions of size I. The dotted straight lines 
show an R^^^ power law, and the dashed horizontal line 
shows the value of the sound speed. The sonic scale can 
be defined as that at which the solid line crosses the sound 
speed level, and is depicted by a vertical dotted line. Be- 
cause of the large error (or, more properly, "scatter") bars 
in the velocity dispersion, the sonic scale is typically de- 
fined only to within factors ~ 4. The caveat should be 
mentioned that the velocity dispersion-size plots we ob- 
tain are clearly not power laws, with progressively smaller 
scales being progressively more deficient in turbulent ki- 
netic energy. This is in part a consequence of the spatially 
varying resolution of SPH, with low resolution at low den- 
sity (in voids) and high resolution at high density (in cores) 
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(Ossenkopf, Klcsscn, & Hcitsch 2001), and in part a conse- 
quence of a power-law range not being expected at scales 
larger than Ad for the small-driving-scale runs. In any 
case, the non-power-law nature of the spectrum at small 
scales implies that our measurements overestimate Ag. We 
discuss the implications of this in §3. 

Figure 2 shows the accretion histories for all runs. It is 
noteworthy that the accreted mass in most runs (except 
the ones with the lowest accretion rates) eventually sat- 
urates at some level. This phenomenon is extremely in- 
teresting in its own right; it clearly appears to be a real 
phenomenon, generally seen in simulations that allow fol- 
lowing the evolution of the system beyond the formation of 
individual collapse sites through the use of sink particles 
(KHMOO; Klessen, Burkert & Bate 2001; Klessen 2001; 
Bate, Bonnell & Bromm 2002). We speculate that it oc- 
curs because the local collapse events decrease the mean 
density of the "interclump" medium, making it more dif- 
ficult to form new collapsing sites from it. Concerns that 
this saturation level be a numerical artifact dependent of 
the resolution are dispelled by fig. 6 of KHMOO, which 
shows no clear trend of the saturation level with resolu- 
tion, but only statistical scatter (different random initial 
conditions with the same global parameters show varia- 
tions of up to a factor of 2 in the final accreted masses). 

Operationally, defining the SFE is not an unambiguous 
task, because not all runs have had time to saturate yet 
(assuming they all eventually do) . Given that recent stud- 
ies suggest that clouds have lifetimes of only a few crossing 
times (Ballesteros-Paredes et al. 1999; Elmegreen 2000), 
we choose to define the SFE as the collapsed mass fraction 
after two turbulent crossing times (tc).^ In the units of the 
code, these correspond to t = 10, 12.5, 6.7 and 4 for runs 
with rms Mach number Mg = 2, 3.2, 6, and 10, respec- 
tively. However, for the runs with Mg = 2 and 3.2, this 
time is longer than they have been evolved. Fortunately, 
the collapsed mass fractions in all these runs have satu- 
rated, and so for them we take the saturation value. Due 
to the statistical variation of the collapsed mass fractions 
mentioned above, we also put error bars of ±50% in the 
reported values of the SFE. 

In order to interpret the results of the simulations, 
fig. 3 shows the SFE versus Aa (a), Mg (6), and Ag (c) 
for all simulations. The linear Pearson correlation coef- 
ficients for the logarithmic data in figs. 3a, 3b and 3c 
are respectively 0.35, —0.73 and 0.75. showing that in- 
deed the SFE correlates best with the sonic scale Ag. 
Note that this is in spite of the fact that the relation- 
ship between the SFE and Ag is clearly not a power law, 
but instead saturates at large Ag. A change of variable 
Ag u = -0.0038 - 0.012 exp -2.39 log Ag gives the fit 
shown by the solid curve in fig. 3c, with a much tighter 
correlation coefficient of 0.86. Nevertheless, the correla- 
tion with Ms is also significant. We discuss this in the 
next section. 

3. DISCUSSION AND COMPARISON WITH PREVIOUS WORK 

Although in the previous section we have fitted a smooth 
curve to the data in fig. 3 c, actually it appears as if there 

are two functional dependences ("modes") of the SFE on 
Ag in this plot: for Ag < 0.03, the SFE is nearly a power-law 

^ We thank the referee for suggseting this criterion. 



on Ag, while it is roughly constant (within a factor of 2) for 
Ag > 0.03. All the points in the former regime correspond 
to cases with Mg = 6 or 10, which are globally supported 
by the turbulence, while most of the points in the latter 
regime correspond to Mg = 2 and 3.2, which are globally 
unsupported. Nevertheless, the large-scale driven (fed = 2) 
high-Mg cases appear to be at the crossover between the 
two regimes. This demonstrates the suggested dual nature 
of the turbulence (Sasao 1973; Falgarone, Puget & Perault 
1992; Vazquez-Semadeni & Passot 1999; KHMOO): glob- 
ally supported runs can be induced to behave as unsup- 
ported ones when the bulk of the turbulent energy is at 
the largest scales. 

It is tempting to speculate that these two "modes" cor- 
respond to the isolated and clustered observed modes of 
SF. In this case, our results suggest that the two would 
not be physically separated modes, but simply represent 
a continuous transition from one functional dependence of 
the SFE on Ag to the other, as the domain transits from 
behaving as if globally supported to behaving as if globally 
unsupported. We see also that, even in the globally un- 
supported case, the efficiency is below unity because of the 
saturation described above. In this case, turbulence would 
only provide the seeds for the internal fragmentation of the 
collapsing molecular cloud region (see, e.g., Klessen et al. 
1998; Klessen 2001), the efficiency being more strongly de- 
pendent on stellar feedback to destroy the cloud. Further 
testing is necessary in order to confirm the correspondence 
between the two regimes found here and the observation- 
ally identified isolated and clustered modes. 

Concerning the resolution of the numerical scheme, as 
mentioned above, Ag is overestimated at small scales by 
our simulations. Moreover, in the absence of numerical 
diffusion, one would expect the density peaks to become 
narrower as well, possibly increasing the number of them 
that go into collapse, and thus increasing the SFE. Thus, 
in general we expect the correlation seen in fig. 3c to be 
shallower if the diffusivity is decreased. Nevertheless, we 
do not expect significant qualitative changes in our results. 

Our reinterpretation in terms of Ag eliminates the ap- 
parent inconsistency implied by previous works that a low 
SFE would seem to require a small driving scale, yet the 
turbulence in molecular clouds appears to be generally 
driven from larger scales. Indeed, our results show that 
a small driving scale alone is not enough to prevent col- 
lapse if it is too weak, while driving at large scales can still 
decrease the efficiency if it is strong enough, as pointed out 
by KHMOO. That is, the controlling parameter is not the 
driving scale nor the rms Mach number alone, but instead 
their combined action expressed as the energy content at 
the sonic scale. Previous studies did not consider large- 
scale-driven cases with a strong enough turbulence level 
to decrease the SFE. 

A few final remarks: 1) The mechanism sketched here 
is expected to operate in addition to any other mechanism 
limiting the SFE based on energy feedback from the mas- 
sive stars, which can disrupt the parent cloud. 2) The oc- 
currence of a saturation of the accreted mass implies that 
clouds could not achieve 100% SFE even if they had arbi- 
trarily large amounts of time available, provided the turbu- 
lence level is kept constant. This is because local collapse 
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uses up part of the mass, rendering the remaining gas more 
diffuse and easy to disperse. Thus, partial dispersion may 
also be part of the natural course of evolution of molecu- 
lar clouds. 3) In this work we have neglected the effect of 
the magnetic fields. Although their effect certainly must 
be incorporated in a full theory, the mechanism discussed 
here is capable of regulating the SFE without relying on 
their presence, depending only on the spatial redistribu- 
tion of the physical quantities in a compressible turbulent 
flow. This is particularly relevant upon recent suggestions 
that molecular clouds may be super- Alfvenic (Padoan & 
Nordlund 1999) and that cores are in general supercritical 



(Nakano 1998; Hartmann, Ballesteros-Paredes & Bergin 
2001; Bourke et al. 2001; Crutcher, Heiles & Troland 
2002). 
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Fig. 1. — Velocity dispersion Av versus size R for all runs considered in this paper (solid line). The nomenclature of the runs indicates 
their nominal rms Mach number Ms and the driving wavenumber fej. The dashed line shows an R}^^ law, the horizontal dotted line shows 
the value of the sound speed, and the vertical line marks the sonic scale. The error bars indicate the scatter of the around the mean value of 
Av. 
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Fig. 2. — Evolution of the collapsed mass fraction for all runs. We define the star formation efficiency (SFE) as the collapsed mass fraction 
at time t = 4, 6.7 for runs with Ms = 10, 6, respectively, and as the saturation value for runs with Mb=3.2 and 2. 



1.00 



a) 



m 0.10 - 

VIIOKB 



0.01 L 
0.1 



Vl:3 2K.1 llgKiXa 

VI6K2 
' :VI10K2 
mK4 



LOG 



0.10 - 



0.01 



1.00 



0.10 r 



c) 



/ 



0.01 



1.0 



10 



0.01 0.10 



Fig. 3. — Dependence of the SFE on the driving length (a), the rms Mach number Ms (b) and the sonic scale Ag in units of the box 
size (c).The error bars in the SFE indicate the statistical variation among different realizations with the same global parameters found by 
KHMOO. 



